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TECHNICAL MEMORANDUM 


COUPLED LOADS ANALYSIS FOR SPACE SHUTTLE PAYLOADS 


I. INTRODUCTION 


This report describes a method for determining the transient response of, and the resultant 
loads in, a system exposed to predicted external forces. For this example, the system consists of 
four racks mounted inside of a space station resource node module (SSRNMO) which is mounted in 
the payload bay of the space shuttle. The applied predicted external forces are Johnson Space Center 
(JSC) forcing functions which envelope worst-case forces applied to the shuttle during lift-off and 
landing. This analysis, called a coupled loads analysis, is used to: 

1. Couple the payload and shuttle models together 

2. Determine the transient response of the system 

3. Recover payload loads, payload accelerations, and payload to shuttle interface forces. 


II. THEORY 

A. Craig-Bampton Reduction 

1. Models 

A finite element model can be thought of as having both internal degrees of freedom (DOF’s) 
and boundary DOF’s where boundary DOF’s connect the model with other models and internal 
DOF’s do not (fig. 1). The models (mass and stiffness matrices) for the racks and the node were 
supplied by Boeing in the Craig-Bampton (CB) format, and the shuttle model was previously built 
and converted to the CB format. The CB format is the result of a CB reduction where the problem is 
reduced in size by ignoring the effects of some of the mode shapes. The size of the problem is 
reduced to the number of boundary DOF’s plus the number of retained mode shapes from an eigen- 
value problem. Because we retain only a portion of the modes for the substructures and for the entire 
shuttle system model, and since the reduced model is in a more diagonalized form, the problem is 
significantly reduced. The following is a discussion of the CB reduction method. 

The motion of the internal DOF’s [X/] can be described as the combination of (1) the internal 
motion with respect to a rigid boundary (all boundary DOF’s are constrained), and (2) the internal 
motion due to the movement of the boundary and the relative deflection of the boundary DOF’s (the 
boundary DOF’s are subjected to different constraints). 


interior dot 


MODEL 

A 



boundary dot 


Figure !. Generic coupled systems. 


2. Constrained Modes 

The first type of motion described above can be analyzed by constraining the boundary 
DOF’s (making a rigid boundary) and then finding the normal modes or “constrained modes” as 
follows: 

1. Constrain the boundary DOF by zeroing out the boundary DOF rows and columns in the 
mass (M) and stiffness (K) matrices. 

2. Solve the eigenvalue problem ( K-w 2 M)PHI n = 0 to get the constrained modes PHl n 
(eigenvectors), and system natural frequencies w (square root of eigenvalues). Only the modes 
below a specified frequency, called a cutoff frequency, are considered significant to the problem and 
are retained in PHI n . Reducing out these unnecessary mode shales reduces the size of the problem 
without losing much accuracy. The determination of the cutoff frequency will be discussed later. 

The internal displacements [X/] with respect to a rigid boundary will now be denoted as 
[Xjcd] for constrain^. We can now write: 

{X kd ) = [PHI n ]{q] , (D 

where q is the modal displacement to be calculated later in the transient analysis. 

3. Constraint Modes 

The second type of motion due to the displacement and deflection of the boundary DOF’s 
(BDOF’s) is analyzed as follows. The equation F = KX can be written in partitioned form such that 
the elements of the matrices which pertain to internal DOF’s can be segregated from the elements 
which pertain to BDOF’s as shown in equation (2), where / denotes internal DOF influence and B 
denotes BDOF influence. Note that Fj = 0 in equation (2) because this equation by definition is 
written to analyze only the effects of movement and relative deflection of the BDOF’s on the internal 
DOF’s (the effect of the internal forces are taken care of in the first type of motion). 
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K„ 

K, b 

1 x , \ 

F, 1 

!=( 0 !. 


k bi 

Xbb . 

i | 

I Fb 

Fb 1 

(2) 


therefore, 

[*//]{*/} = 0 
{X/} = -[K u ]- l [K m ] {X B } , 

if we let 

-[K„r x [K IB } = [PHId 


then 

{X;} = [PHI C ] {X B } . (3) 

To gain more insight into what PHI C is, we can let X B be an identity matrix. In this case, each 
column of PHIc is a map of relative internal displacements (mode shape) resulting from a unit dis- 
placement of one BDOF while fixing (constraining) all the other BDOF’s. For this reason, we call 
PHIc the constrain! modes, and we will now refer to this type of X/ as X/ c ,. Rewriting equation (3) 
with this new notation we get: 


[X lcl ] = [PHI C ] {X B } . (4 

Xic is the matrix of internal displacements due to the relative displacements of the BDOF’s found in 
X B . 

4. Equations of Motion 

The complete internal motion is found by combining the constrained and the constraint types 
of motion. 


{X/} = { X lcd ) + {X Ict } . 

Substituting equations (1) and (4) into equation (5) we get: 

{X/} = \PHI n ]{q}+[PHld{X B } 


and since 


{X B ) = {X B } , 


we can write: 


X, 

X B 


PHI n PHIc 
0 1 


(5) 

( 6 ) 


(7 1 


q 

x B 


( 8 , 


We now define 


PHI n PHlc 
0 1 


= [TCB] , 


( 9 ) 


where [TCB] is the CB transformation matrix. 

The number of columns in PHI ^ is equal to the number of retained modes, and the number of 
columns in PHI C is equal to the number of BDOF’s. Therefore, the size of [TCB] is equal to the 
number of retained modes plus the number of BDOF’s. We will later see that the size of our reduced 
model is this same size. 


The equation of motion for a system with forced vibration and no damping (damping will be 
taken care of later) is: 

[M]X + [K]X = {F} . (10 ) 

If we write this equation in partitioned form with respect to internal and boundary DOF’s, we can 
get: 

Mil M ib I Xi 

. Mbi Mbb J \ x B 

where I - interior DOF’s and B = boundary DOF’s. Mu is a matrix of those mass matrix terms 
which represent only internal DOF’s. Mbb is a matrix of those mass matrix terms which represent 
only BDOF’s. M jb and Mbi have terms with a mixture of boundary and internal DOF influence. Ku, 
K/b, Kbi, and Kbb are defined the same. By substituting equations (8) and (9) into equation (11) we 
get. 


1* 


K„ K lB 
[ Kbi Xbb 




( 11 ) 


Mu M/b 

Mbi Mbb J 


[TCB] 


I * 

\ X B 


K„ K,b 
Kbi Kbb 


[TCB] 


<1 U I Fi 


X B 


F b 


(12) 


To simplify and to make symmetric, we multiply through by [TCB] T and get, after several 
steps and substitutions. 


/ W es ] I « 1 + 

M B q Mbb J \ x B j 


w 2 0 

0 Kbb J 


Q 

X B 


= [TCB] 1 


tI Fi 

Fb 


(13) 


where 


Mqb = [PHI N ] T [M„][PHI C ] + [PHI N ] T [M, B ] 
M B q = [PHI c ] T [M„][PHI N ]+[M B i][PHI N ] 


w 2 = [PHInY [K//][PHIh] (squared constrained system natural frequencies) 
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If we let 


1 m QB = [MCB] 
Mbq Mbb J 


and 

w 2 o 

= [KCB] , 

.0 K B b. 

then we can rewrite equation (13) to get the CB equation of motion: 

[MCB] [ 4 UlKCfi]! « j = [TCB[ t \ p j. 

\ iig I I x ® I Fb I (14, 


B. Coupling 

To couple systems together, we put the models (mass and stiffness matrices) into a common 
coordinate system and then add the equations of motion. 

As an example (which directly applies to the node problem), we will couple two subsystems 
called rack (R) and node (AO into a new system called the payload (P) system. Then we will couple 
this payload to the shuttle (5) system to get the complete lift-off ( LO ) system. The R system 
attaches to the N system, and the P system attaches to the 5 system at the node (fig. 2). 

When CB reductions are done, the BDOF’s must include all boundaries to be used in future 
couplings. For example, the node attaches to the rack and to the shuttle. Therefore, the BDOF’s for 
the node include the node to rack ( n/r ) and node to shuttle (n/s) BDOF’s. 

Since the models of the R, N, and S systems are supplied in the CB format, we are given the 
CB transformations [TCBr], [ TCB #], and [TCBs], where 

j XlN \ j QN \ 

(*/*) = [TCB R ] ! CfR I , and ' Xb n,r { ~ WCBn\ ( Xb n/r / > ( 15 ) 

\Xb r/n | I Xb r/n I * Xg n/s I ' Xg n/s I 

and 

j Xis J = [TCBs] j qs I 

l Afl sin j | A# s/n I 

Note that n/s stands for node to shuttle coordinates and s/n stands for shuttle to node coordinates. 
These TCB's, which convert from modal coordinates to discrete coordinates, were generated as 
shown above during the initial R, N, and S CB reductions. To couple the R system to the N system, 
we first write the equations of motion from equation (14) for the rack and then the node systems. 
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Figure 2. Lift-off coupled system. 


[MCB] r ( qR | + [KCB] r J qR j = [TCB] t r (/'*), 

\X B J \ FBr ' n ' 


(16a) 


| 4 n 1 

1 qN 

T 1 

( F,n \ 

XBn/r 

1 + [XCB] N ( X B n/r 

: = [TCB) r h 

j ^ B n/r I 

Xb n/s 

\ X B n/s 

1 1 

l Fb n/s 1 


(16b) 


To couple the rack and node systems together, we need to convert them to a common coordinate 
system which contains all coordinates found in the rack and node systems. Such a system, in this 

I <1r 


case, is the 


qN 

X B r/n 

X-B n/r 
Xb n/s I 


system. Since we are coupling the rack and node at the r/n boundary, we are 


qR 


imposing the relation that Xb r/n = Xb n/r. Therefore, we will convert to a j q n \ system, which 

Xb rin I 


X B 


nfs 


also contains all coordinates found in the rack and node systems. 
We can write: 


qR U 

Xb r/n 


= [TR] 


I I 

I qN \ 


qR 

qN 

Xb r/n 
Xb n/s 


qN 

, and { X Br/n 

Xb n/s i 


qR 


= [TN]{ & 

^B r/n 
Xb n/s 


(17) 
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where [7/?] and [77V] are the coordinate transformation matrices which convert from the new com- 
mon (coupled) coordinate system to the rack and node systems, respectively. Unless a coordinate 
rotation is involved, these transformations are generally made up of 0 s and 1 s. From equations 
(17) we can solve for [7/?] and [77V] which we will need later in this analysis to convert back to the 
component system coordinates (R and TV) from the common coupled system coordinates when 
recovering component loads. In this example. 


[77?] = 





' 0 

1 

0 

o' 

10 0 

0 ' 

and [77V] = 

0 

0 

1 

0 

.001 

0. 

. 0 

0 

0 

1 - 


Substituting equation (17) into equations (16a) and (16b), we get equations i.18) and (19). 


[MCBMTR] 


QR 

'4n 

Xb r/n 
Xb n/s 


+ [ KCBMTR ] 


<7* 

qN 

Xb r/n 

Xb n/s 


=[tcb] t r \*: ir , , 

TB r/n 


(18) 


[MCB] n [TN] 


<IR 

QN 

Xb r/n 

Xb n/s 


+ [ KCB] n [TN] 


( qR 1 

1 

f fin \ 

1 qN ' 

1 X B r/n 
\ Xb n/s i 

; = [TCBf N j 

j F B n/r ) 

( Fb n/s i 


(19) 


Now we will multiply equations (18) and (19) by [TR]T and [TN] t , respectively, to simplify 

I 


[TR] 1 [MCBMTR] 


<1R 

qN 

Xb r/n 
Xb n/s 


QR 


+ [TR] t [KCB] r [TR] l qN )= [TR] t [TCB] t r 

*B r/n 


\ 


r_ I Fir 

Fb r/n 


, ( 20 ) 


Xb n/s 


[TN] t [MCB] n [TN] „ 


qR 

qN 

Xb r/n 
Xb n/s 


qR 

q A 

1 B ‘ 

\ Xb n/s 


+ [TN] t [KCB]n[TN] qN = [TN] t [TCB] t n Fsn/r (21) 

1 Xb r/n 


Fin 

Fb n/r 
Fb n/s 


Now, we can add equations (20) and (21) while letting 

[TR] t [MCB]r[TR] + [TN] t [MCB] n [TN] = [MP] , 


(22a) 
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and 


[TR] t [KCB] r [TR] + [ TN] t [KCB] n [TN] = [ KP ] , 


(22b) 


where MP and KP are the P mass and stiffness matrices. The result is the general coupled payload 
equation of motion: 


[MP] { 


<?* 


1 <7* \ 

4n 

! + [KP] ! 

<7 N 1 

Xg r/n 

1 Xg r / n j 

1 ' 

l Xg n/ s ) 

Xb n/s , 




= [TR] t [TCB] t r j + [TN] t [TCB] t n | Fair J • (23) 


It is important to note that for shuttle payloads, there are generally no internal forces applied 
to the rack or the node, so F IR = F ]N = 0. After making the substitutions for [77?] r , [TCB] t r , [TN] t , 
and [TCB] t n , the right hand side of equation (23) simplifies as follows: 


1 O' 

r 'r 

. 



'0 0 o' 

r 



0 0 

phi t n 

0 

( 0 

+ 

1 0 0 

PHI l N 

0 1 

0 1 

L PHI T c 1r . 

\F B r/n 


0 1 0 

PHI T c 

In ' 

0 0. 





.001- 





I 

1 0 

| 

f 

0 



0 \ 


J 

0 

u 

1 

0 


^ 

0 



[F B n/r)\ 

[Fb r/n) + { Fsn/r } j 


0 


1 

[ { nls] i 1 

1 

{ F b n/s 

) 


i { Fb nfs } | 


0 

F B n/r 
Fb n/s 


L 


0 
0 

{ Fb rln } 
0 


The last step in the simplification is due to the fact that the sum of the forces at the coupled boundary 
between the rack and the node are zero, so Fg „/ r +Fg r / n = 0. This assumes that there are no 
externally applied forces at the boundary. Using the above simplifications, equation (23) becomes 
equation (24), the final payload equation of motion. 


[MP] 


I 4R 


' <IR 


0 1 

'4n 

+ [KP] i 

qN 


0 1 

Xb r/n 

I Xg r / n 


0 

i Xg n/s 1 


l Xg n / s 


, Fgn/s ) 


(24) 


To further illustrate that this method couples by simply adding at common BDOF’s, we can 
expand equation (22a) to get equation (25) and then multiply to get equation (26). 


10' 


0 0 

/ Mqb 

0 1 

_MgQ Mbb Jp- 

.0 0. 

A 


0 

0 


0 

1 


0 

0J 


0 

1 

0 

0 


0 

0 

1 

0 


o' 



0 

I Mqb 


0 

Mbq Mgg 

N _ 

1 . 



1 

0 

0 


0 

1 

0 


0 

0 

1 . 


= [MP] 


( 25 ) 
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= [ MP ] . 


(26 


I 0 MqBt 


r 0 

0 

0 


I 0 

MqBt 

0 0 0 

+ 

0 

r / 

Mqb^\ 

— 

0 I 

Mqbh 

_Mbqt 0 Mbbt. 


0 

l M B Qn 

M B Bn J 

J 

M B Qr M B Qn 

(MsBr+MaBn) 


Expanding equation (22b) would show the stiffness matrices coupling in the same way. 

We will next couple the payload to the shuttle. To simplify the problem, we will do a CB 
reduction as before on the new payload system, saving only die boundaries needed for coupling; 
those common to the payload and shuttle systems. The CB reduction, performed as before, gives us 
a new CB transformation matrix [TCB] P which satisfies the following relations. 


[TCB] P = 





QR | 



phi n 

PH1 C ' 

, and 

<7;V ' 

} = [TCB]p 

f qp 

0 

I 

P 

X B r/n 
, X B n/s 

1 Xb n/s 


(27) 


Again, we can get the CB mass and stiffness matrices. 

[MCB] p = [TCB] t p [MP][TCB] p , (28a) 

[KCB]p = [TCB] T p[KP][TCB] P . (28b) 


The new CB equation of motion for the payload is: 

v qp \=[TCB] T p[TN] i 

*Bn/s ) 


[MCB] P j qP \ + \KCB)p 

1 X Bn , s I 


0 

Fb n/s 


(29) 


The CB equation of motion for the shuttle is derived the same way and is: 


[MCB] S 


qs 

Xb s/n 


| + [KCB]s j ^ | = [TCB] T s 


Fis 

F B s/n 


(3(0 


Since there are externally applied forces on the shuttle internal DOF’s, the F/ term does not go to 
zero as for the payload system. 

To couple the P and the S systems into the LO system, we must rearrange and add as we did 
for the rack and node subsystems. We can derive the coordinate transformation matrices [TP] for 
the payload, and [FS] for the shuttle from equations (31). The LO mass and stiffness matrices 
([ MLO ] and [ KLO ]) can then be found by using equations (32a) and (32b). 


( 

1 4P 


j qs 
| X B n/s 

\ 

qp | 

qp 

\ X B n/s 

-H • i 

1 \ n/s 1 

■ 

= [TS ] | 

qs | 

l Xb n/s I 
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[TP] t [MCB] P [TP] + [TS] t [MCB] s [TS] = [MLO] , 


(32a) 


and 

[TP] t [KCB] p [TP] + [TS] t [KCB\s [TS] = [KLO] . (32b] 

As before for the payload system, we substitute equations (31) into equations (29) and (30). We 
multiply these modified versions of equations (29) and (30) through by [TP] T and [TS] T , 
respectively, and then add them together. Substituting equations (32) into this new equation we get 
a coupled LO equation of motion. 


[MLO] 



= [ TSY[TCB] T s 


F/s 

Fb sip 


[TPflTCBfpiTN] 1 


0 

Fb n/s 


(33) 


Note that the payload attaches to the shuttle at the node (fig. 2). This means that F B „/ s is the same 
as F Bp / s . Since the payload and shuttle are coupled at the common boundary, then F Bs/p + F Bp / s = 0 
(again assuming no externally applied boundary forces), and therefore F B s/ P +F B nJs = 0. Using this 
relationship after substituting all of the transformation matrices into equation (33) and then simplify- 
ing, we get equation (34), the complete LO equation of motion. 


<iP 


qp 


[MLO] [ q s J + [KLO] | qs 

X B n/s I 


X B 


n/s 


= [TS] t [TCB] t s \ Fi o s 


(34) 


This is the basic equation used in the transient response analysis (see below) to get the modal 
responses (qs) from the forcing function (F/ s ). 

C. Transient Response 

Equation (35), the eigenvalue problem, is solved for the eigenvalues ( w^o 2 where wpo are 
the system free-free natural frequencies), and the eigenvectors ( PHI LO ) which are the system mode 
shapes. 


where 


([ X lo ]-[w lo 2 ] [ M lo ]) [PH1 lo ] = 0 , 


I qr I 

qs ) = \PHI LO \[qio) ■ 

' Xfi n/s I 


(35) 


(36) 


Since the mode shapes ([PHI LO ]) are just relative amplitudes, they can be multiplied by a constant 
and still retain the same shape. We can normalize the mode shapes by picking constants such that 
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[PHI L of[M LO ] l PHI lo ] = [/] • 


(37) 


When equation (37) is satisfied, then so is equation (38). 

[PHI lo ] T [K lo ] [PHI lo ] = [w lo 2 ] • (38) 

Now equation (36) can be substituted into equation (34) and, after multiplying through by 
[PHI lo ] t to simplify, we can use equations (37) and (38) to get the normalized equation of motion 
with no damping. 

[/] (wol + I'Vz.o 2 ] Ulo) = [PH1 lo ] t [TS\ t \TCB] t s j j • < 39 ' 

At this time it is convenient to take damping into account. Damping is assumed to be propor 
tional damping [2] as defined by equation (40). 

[D] = [2 zw) , (40) 

where z is the damping ratio, which is assumed for lift-off to be 1 percent of critical damping for 
modes below 10 Hz and 2 percent of critical damping for modes above 10 Hz. For landing analysis, z 
is assumed to be 1 percent of critical damping for all modes. Since damping force is proportional to 
velocity, our final equation of motion is then: 

[I] {qLo} + [2z w Lo]{qbo} + [ w L0 2 ]{qLo} = [PHIlo\ T [TS] t [TCB] T s^ j • (41) 

Note that we now know everything in this equation except the q’s and the external forces on 
the shuttle applied to the internal DOF’s ({F/s}). The forces are supplied by JSC in the form of 
forcing functions which are predictions of a certain shuttle environment during lift-off or landing. A 
forcing function is a matrix where each column represents specific forces on specific shuttle DOF’s a 
one instant of time. The number of columns then is the number of time increments to be studied. Lift 
off is usually run for about 10.5 s, and landing for about 2 s. The time step used can vary, but 0.005 s 
is typical. 

By virtue of the modal analysis, equation (41) is a set of n linearly independent differential 
equations with n unknowns ( n = size of the model). We can now solve for the q s for the first column 
(time step) of the forcing function. The results of this solution are used as the initial conditions for 
the solution of the next time step. The solutions keep “feeding" each other like this until all forcing 
function columns are used. This gives us the values of q for all time steps during an event such as 
lift-off. The q s represent the modal displacement, velocity, and acceleration of the system. 

D. Recovering System Response From Modal Response 

Now that we have solved equation (41) for the modal displacement and acceleration 
responses of the system ( qs ), we need to determine what this means with respect to specific inter- 
face forces between subsystems and net center-of-gravity (CG) accelerations of subsystems. The 
hybrid modal acceleration method is frequently used to get interface forces and CG accelerations 
when the system is either determinate (six BDOF’s) or slightly indeterminate (seven BDOF’s). In 
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some instances, even seven BDOF’s require another method called the acceleration method, which 
is outside the scope of this report. The following makes use of the hybrid modal acceleration method. 


The rigid body displacements (first six rows of q) are set equal to ziero which forces the 
coordinate reference point to be the CG of the shuttle. Otherwise, the large movement of the entire 
shuttle system would be included in the displacement results. The modal displacements and accel- 
erations are then transformed back to the point of interest, as shown below, using the transformation 
equations derived above and repeated here for clarity. 


and 


I qr \ 

\ IS = [PH1 W ] . 

\ Xb n/s I 


qp 

Xg n/s 


I qp \ 
[7V>1 « • 

\ *B n/s I 


qs 

Xb n/s 


I qp \ 

= [TS] qs 

1 Xg n/s I 


[TCB] P = 


PHI n PH1 c 
0 / 


qR \ 

Xb r/n I 


= [77?] 


qR 

qN 

Xb r/n 

Xb n/s 


and 


\ qR | [ 

, and 


qR 

qN 

Xb r/n 

Xb n/s 


qN 

j ’ al,u \ Xb r/n I 

I l X s n/s I 


= [TCB]p 


\ = [ TN ] 


qp 

Xb n/s 


qp 

qN 

Xb r/n 

Xb n/s 


\ 

r 


X,R 

Xb r/n 


= [TCB R ] M* 


a: 


B r/n I 


and 


X/n 

Xgn/r 

Xb n/s 


I qN 1 


qN 

= [TCB N ] ; x Bn /r I 

1 X B n/s I 


(36) 

(31) 

(27) 

(17) 

(15) 


j X ls 

1 Xb s/n 


[TCB S ] 


qs 

Xb s/n 


These equations can be used to find the response of any of the subsystems. For example, if 
we want the response of the rack during lift-off, we back into it from the system modal response as 
follows: 


I qp 

\ » 

' Xg n/s 


= 1 PHhxtMqijo) . 


qp 

Xb n/s 


= [TP] 


qp 

qs 

Xb n/s 


\ 

r 
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q R 

W } = [ TCBp ] 

J Xb nls I 
\ Xb r/n I 


4P 

Xb n/s 


I q R 

r/n 


= [77?] 


q R 

q N 

Xb r/n 
Xb n/s 


j XI R 

r/n 


= [ tcb r ] 


q R 

Xb r/n 


Therefore, the modal response of the rack is: 

j « j = [TR][TCBp][TP][PHI LO ]{qLo) , 

\ Xb r/n | 


and the rack discrete response is: 

j XlR \=[TCB R ][TR][TCBp][TP][PHI LO ]{qw ) 

l-^fl r/n 

Similarly, the modal response for the node is: 


I L 

1 y" r 

\ n/s I 


[TN][TCBp\\TP]\PHIlo\ {quo) . 


and the node discrete response is: 


( X IN ) 

X Bn / r ) = [TCB N }[TN][TCB P ][TP}[PHI LO ]{qLO ) 
\ X Bn/s I 


E. Recovering Loads Using Loads Transformation Matrices 

Now that we have the response, we need to convert this information into loads 
use of loads transformation matrices (LTM’s). The LTM’s for finding interface forces 
center-of-gravity accelerations (CGA’s) are derived in the following. 


(42a) 


(42b) 


(43a) 


(43b) 


through the 
(IFF’s) and 
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1. Interface Forces 


(13): 


IFF’s are simply the boundary forces (F B ) between two subsystems. Recalling equation 


I M 


QB 


M B q Mbb] 


I ? 

\ X B 


w 


0 


1 1 ■ |rCB| i Z 


0 Kbb j 

There are no applied internal forces in the rack or node (F f = 0) and since 

[TCBf = 

equation (13) simplifies to: 


PHI n 

PHI C 


0 

I 


I M, 


QB 


M B q Mbb 


I <? 

\ x B 


w 


0 


0 K bb 


Q 

X B 


from which we see that the boundary forces (same as IFF’s) are: 


0 

Fb 


Liff = [ Fb ] = [M B q M B b] { } + [0 Kbb] ( J 

I X B I \ X b 


( 13 ) 


(44) 


(45) 


[M B q M B b] is sometimes called LTMA, the acceleration LTM, and [0 Kbb] is called LTMD, the 
displacement LTM. Rewriting equation (45) we get equation (46): 


L,ff = [LTMA] 


I i 1 

{ * J + [LTMD] 
\ X B I 



(46) 


Substituting equation (42a) from the example above into equation (46), we get equation (47) which 
gives the interface forces between the rack and the node: 


Liff r = [ LTMA] r [TR][TCB p ][TP][PHI lo ] [q LO ] 
+ [LTMD] R [TR][TCBp][TP][PHI LO ] { qw } - 


2. Net CGA’s 

A net CGA is the average acceleration of a subsystem CG assuming the subsystem to be a 
rigid body. This should not be confused with the acceleration of a point at the CG. We need to know 
the forces at the CG of the subsystems. Since we now know F B , the forces at the boundary 

(equation (45)), we need to transform them to the CG. To do this we solve for the rigid body trans- 
formation (RBT) such that: 


[X B ] =[RBT] {X CG } 


(48) 
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( 49 ) 


Since [RBT\ is orthogonal, then [RBT] T = [RBT] l , and it follows that 

{X CG ) = [RBT\ t {X B ) , 


and that 


l Fee ] = [RBT] t [F b ] 


(50) 


To convert equation (45) to the CG coordinate system, we multiply through by [RBT] T and get: 


[RBT] t [F b ] = [RBT] t [M bq M bb ] | * | + [RBT) t [0 K bb ] 


q 

x B 


(5i) 


Since converting to the CG coordinate system is an RBT, all coordinates retain the same relative 
position. Therefore, for CGA’s we assume a rigid model which means there are no forces due to 
deflections and so: 

[RBT\ t [0 K bb ] = 0 . (52) 

Substituting equations (50) and (52) into equation (51) we get: 


[F C g) = [RBT] t [M bq M bb ] 


X B 


Since 


then 


[ Fcg] = WcG ] {Xcg) > 


{ X CG } = [McgT 1 [Fcg] ■ 

Substituting equation (53) into equation (55) we get: 

[X CG ) = WcgY 1 [RBT] t [M bq M bb ] { * 


(5:5) 


(54) 


(55) 


X B ) 

We can get [Mcg\ by performing the transformation: 

[M C g] = [RBT] t [M bb ] [RBT] . 
Substituting equation (57) into equation (56) we get: 

[Xcg) =([RBT\ t [M bb ][RBT\T 1 [RBT] t [M bq M bb \[ ? 

I X B 


(56) 


(57) 


( 58 ) 


If we let 


[ACG] = ([RBT] T [M BB ][RBT)r 1 [RBT] T [M BQ M BB ] , (59) 

then from equation (58), 

{Xcg)=[ACG]1 5 1. 

' * B ’ (60) 

Substituting equation (42a), from the example above, into equation (60) we get equation (61) which 
gives the CGA’s for the rack during lift-off. 


{X C g)r = [ACG] r [TR][TCB p ][TP][PHI lo ] {quo) ■ 


(61) 


III. SUMMARY 


Coupled loads analyses require a lot of bookkeeping techniques. First the models are coupled 
while creating a trail of transformation equations in order to later back out to the subsystem level. 
The equations of motion are then written for the coupled system. Next, the transient analysis is 
performed by applying the forcing functions to the equations of motion to get the system displace- 
ments and accelerations during the dynamic event. Finally, the system transient results are used 
along with the transformation equations to back into specific subsystem responses. At this time, it is 
important to point out that we can calculate many more different kinds of component responses other 
than interface forces and CGA’s. These responses were needed for the SSRNMO and therefore 
were used in this report. The topic of a future report might be the usage of LTM’s to recover other 
types of responses. 
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